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Many-body stochastic processes with weighted multiplicative interactions are investigated ana- 
lytically and numerically. An interaction rate between particles with quantities x,y is controlled 
by a homogeneous symmetric kernel K{x,y) oc x'^y^ with a weight parameter w. When w < 0, 
a method of moment inequalities is used to derive log-normal type tails in probability distribution 
functions. The variance of log-normal distributions is expressed in terms of the weight w and inter- 
action parameters. When interactions are weak and a growth rate of systems is small, in particular, 
the variance is in proportion to the growth rate. This behavior is totally different from that of 
one-body stochastic processes, where the variance is independent of the growth rate. At w > 0, 
Monte Carlo simulations show that the processes end up with a winner-take-all state. 
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A number of probability distribution functions (PDF) 
with fat tails of log-normal type have been observed 
in a wide variety of processes including not only nat- 
ural but also social and economic phenomena^, 0, 
When we conduct a keyword search on " log- normal" , we 
find thousands of literatures in various fields, say, size 
distributions of small grains 0, |^ 0, H, aerosols 
cloudsUl m [13, foamsPj, galaxieslll E [13, the 
abundance of species 0, JwL 01, the number of em- 
ployees in manufacturing plants|2^, the prices of in- 
surance claims [21I li^ . and the size of farms [23, and 
so on. In spite of accumulations of a vast amount 
of data, understandings of log-normal distributions are 
still quite poor compared with those of power-law dis- 
tributions. For example, the concepts of self-organized 
criticalitylH, [H, [13 and scale-free network[23 113 are 
useful in explaining power-law distributions. For log- 
normal distributions, in contrast, there exist no corre- 
sponding explanations. It is well-known that in systems 
of a single degree of freedom, a multiplicative stochas- 
tic process exhibits a log-normal distribution[13, [13- In 
most of actual processes, however, interactions play a sig- 
nificant role and systems must be treated as those with 
multiple degrees of freedom. In addition, experiments 
and simulations give only circumstantial evidences. So, 
there is a real need for analytical derivation of log-normal 
distributions in many-body systems. To our knowl- 
edge, however, the theory of log-normal distributions has 
not been developed successfully in many-body processes 
which is irreducible to one-body problem. 

In recent years, stochastic processes with multi- 
plicative interactions have been attracted considerable 
attention[33, [^ [13, [13 ■ Although the processes are 
simple extensions of the one-body multiplicative stochas- 
tic process to many-body processes, it does not exhibit 
a log-normal distribution but a power law at the tail 
of PDf[13. In this Letter, we investigate many-body 



stochastic processes with weighted multiplicative inter- 
actions analytically. In consequence, we succeed in de- 
riving a log-normal type tail. Furthermore, we find that 
when interactions are weak and a growth rate of sys- 
tems is small, the variance of log-normal distribution is 
in proportion to the growth rate. This behavior is to- 
tally different from that of one-body systems, where the 
variance is independent of the growth rate. 

We consider a system of N particles with positive 
quantities Xi{> 0) (i = 1, . . . , N). At each time step, the 
system evolves with a binary interaction between parti- 
cles labeled by i and j {i 7^ j), and two quantities Xi, Xj 
are transformed into a;^, x' by the rule 
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I3xj 



(3xi 
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where a, (3 (> 0) are positive interaction parameters. In 
the limit N ^ 00, the processes are described by the 
master equation 

df{z,t) _ 



dt 



dx / dy f{x,t) f{y,t) K{x,y;t) 
Jq 

X i [ d{z — {ax + f3y)) + S{z — (J3x + ay)) 
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where f{x,t) is a PDF of the quantities and K{x,y) is 
a kernel representing interaction rates. It has been al- 
ready reported that when a kernel K is constant, in 
other words, the interaction takes places between ran- 
domly chosen particles, the PDF has a power-law tail[13- 
The exponent of the tail is a continuous function of pa- 
rameters a, (3 and is calculated analytically via a tran- 
scendental equation. In many cases, however, the kernel 
K depends on quantities x and y. A typical example 
is the case where a quantity x represents a mass M or 
size R of particles. Generally, a transport (diffusion) co- 
efficient and then an interaction(reaction) rate depend 
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on M and R. So, it is natural to think what happens 
in the processes with non-constant kernels. To answer 
this question, we address the case with a homogeneous 
symmetric kernel 



K{x,y\t) 



(3) 



where u; is a weight parameter and rriw is a w— th order 
moment defined by mw(t) = f{x,t)dx. Hereafter, 

we deal with the case a > 1 or /3 > 1. In this case, the to- 
tal sum of the quantities increases and the system grows 
as a whole. When w is positive, interactions between par- 
ticles with large quantities are accelerated. Thus, the dis- 
tribution is considered to become wider than the power- 
law. At w < 0, the PDF is expected to have a narrower 
tail. Firstly, the case w < is investigated by considering 
an asymptotic moment function of PDF. A log-normal 
type tail in PDF can be derived analytically and numer- 
ical results support this finding. Secondly, Monte Carlo 
simulations are performed in the case w > 0. The exis- 
tence of winner-take-all type distribution is elucidated. 

As mentioned above, we treat the case where the sys- 
tem grows and does not have a steady state. Then, we 
attempt to find a scaling solution of Eq. Q and assume 
the relations 



7* 



(4) 



where 7 is a scaling parameter representing a growth rate 
of systems. Substituting Eq. |@J into Eq. we have a 
scaled form of the master equation. 
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where /ip is a scaled th moment. The parameter 7 is 
defined by Eq. © at p = 1 , 



7=(a + /3-l) 



(6) 



Note that Eq. (0) is symmetric with respect to interaction 
parameters a and /3. Without loss of generality, there- 
fore, we put a > [3. To proceed further, we follow the 
method of moment inequalities which is used to study 
inelastic hard sphere models by Bobylev et a/.js^]. Ap- 
plying lemma 2 of Is^l to Eq. lO, we obtain inequalities 




where kp denotes the integer part of (p -I- l)/2 for p > 1. 
It should be noted that when p is an integer n e N, 
inequalities Q) are reduced to the equality 



n{a + P-l)fin = (a"+/3"-l) 



(8) 
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The goal of this analysis is to estimate the asymptotic 
form of moment function /Up (p —^ 00) when w < 0. Here, 
we assume the asymptotic form of p— th moment 



Hp — exp(ap^ — bplnp + cp), 



(9) 



where a, b, and c are unknown constants. It is found that 
Eq. lO satisfies inequalities 



0< ( n Mfc+^Mp^fc+.. <Cexp{-Ap + B\np] ,{1Q) 



for 1 < A; < fcp, where A,B,C{> 0) are certain positive 
constants. It follows that in Eq. (0, l.h.s. of the first in- 
equality and r.h.s. of the second inequality vanish in the 
limit j3 — !■ cx). Hence, the following recurrence equation 
is fulfilled in the p— asymptotic region 



p{a + l3- l)^lp = (aP + 13^-1) 



-, (p^oo).(ll) 



Substitution of Eq. 10 into Eq. ((TT|l leads to 
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where 5ap is the Kronecker delta. Consequently, the 
asymptotic solution of moment functions are determined 
self-consistently. We compare this result with data of 
Monte Carlo(MC) simulations and numerical solutions of 
the recurrence equations © and Hll|) . Values obtained 
by MC simulations are used for the initial values of the 
recurrence equations. Notice that Eq. ^ and pro- 
vide exact results asymptotically. Two typical results at 
w — —1 and —0.5 are illustrated in Figs.^ In the large p 
region, ln(/ip) behaves as a quadratic function with a cur- 
vature consistent with a in Eq. 1)12(1 . MC results agree 
well at small p. Deviation at large p comes from the 
finiteness of N . 

Next, let us consider the shape of PDF. Inverse Mellin 
transform of the asymptotic moment function Eq. 
with lO 



1 l'l-\-ioo 
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FIG. 1: Semi-logarithmic plots of the scaled p— th moment fip 
versus at a weight parameter u; = — 1 (a) and —0.5 (b). 
Interaction paramaters are a — 1.43589 and /? = 0.1. Data 
of MC simulations are illustrated in solid lines. In MC simu- 
lations, the number of particles is A'' = 10^ and the number 
of interactions is T = 200 x A''. Solutions of the recurrence 
Eq. JHJ and lllll are plotted in square dots ((a) only) and 
triangular dots, respectively. The results of Eq. Q with pa- 
rameters 1121 are shown in dashed lines. 
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FIG. 2: (a) Semi-logarithmic plot of temporal evolution of the 
first order moment mi{t). The inset shows the divergence of 
mi(t). (b) Double logarithmic plot of temporal evolution of 
the PDF. {a = 1.09987492, = 0.005, and w = 0.5.) In MC 
simulations, the number of particles is A'^ = 10^ and unit of 
time is defined by the number of interaction steps A'^/2. 



gives PDF. It is immediately found that the leading order 
term of Eq. (jSJ is transformed into a log-normal distri- 
bution. 



■ exp 
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It is also checked by MC simulations that the PDF is log- 
normally distributed. It should be emphasized that the 
variance of the log-normal distribution is determined by 
a — ln(a)/2|ii;| and is independent of the smaller inter- 
action parameter (3 (< a). In many actual phenomena, 
in addition, interactions are weak and a growth rate is 
small. In this case(0 < ^9 < 7 < 1), Eq. I© and lO 
give 



1 ^liH^u 



-7 cx 7. 



(15) 



2\w\ ^ii+w 

This dependence of the variance a on the growth rate 



7 differs qualitatively from that of one-body processes, 
where a is independent of 7. This shows that the log- 
normal distributions of many-body processes belong to a 
totally different universality class from those of one-body 
processes. 

The sign of the weight parameter w is essential to the 
tail of PDF. When = 0, Eq. (|ll|l reduces to a tran- 
scendental equation (a -|- /3 — l)p = -f /J^ — 1, which 
coincides with the result in [s^l- When w > 0, it is sug- 
gested by MC simulations that any scaling solution does 
not hold. This is confirmed by the fact that the growth 
rate 7 is not kept constant and mi(t) diverges drastically 
after a critical time tc as shown in Fig.[5fa). Instead, in 
the long time limit, the processes end up with a winner- 
take-all state where two particular particles are selected 
almost invariably and gain almost all of the total sum 
of quantities as illustrated in Fig. |2tb) . It becomes evi- 
dent that the power-law distribution(i/7 = 0) emerges at 
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the boundary between the log-normal tail(w < 0) and 
the winner-take-all state(w > 0). This situation is sim- 
ilar to those of self-organized criticality and scale-free 
network [23. 125I EgI . that is, criticality appears at 

the edge of chaos |35ll36 

As far as we know, this work gives the first analyti- 
cal derivation of log-normal type distribution in many- 
body processes which is irreducible to one-body prob- 
lem. It can be concluded that the following two condi- 
tions are essential for the emergence of log-normal dis- 
tributions: (i) systems grow with multiplicative inter- 
actions, (ii) interactions between particles with smaller 
quantities are accelerated. In this Letter, we have treated 
only the case with a kernel given by Eq. (jSJ. Recently, 
we have examined the processes with a generic kernel 
K{x, y) oc x~'^y~'' + x~^y~"' (a, 6 > 0) and obtained es- 
sentially the same results. The details will be reported 
elsewhere jsT]. 

The processes are rather simple and expected to be- 
come a prototype of certain log-normal type distribu- 
tions in nature. For example, production methods of fine 
particles are categorized into two types: breaking-down 
and building-up. The former may be described by one- 
body processes as the zeroth approximation. However, 
the latter must be treated as many-body processes be- 
cause interactions between particles are essential for the 
building-up process in general. This suggests that log- 
normal distributions by building-up methods belong to 
the universality class of many-body processes with a oc 7, 
while those by breaking-down methods are in the class of 
one-body processes with 0^7'' = const. The difference 
is expected to be verifiable experimentally. 
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